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In inhomogeneous bulk ferromagnets, the dominating sources of spin disorder are related to spatial 
variations of (i) the magnitude of the local saturation magnetization and of (ii) the magnitude and/or 
direction of the magnetic anisotropy field. For the particular example of a porous ferromagnet, 
where the magnetization inhomogeneity is at maximum, we demonstrate, by means of experimental 
neutron scattering data and micromagnetic simulations, the anisotropic character of magnetization 
fluctuations induced by the dipolar interaction. 


I. INTRODUCTION 

In polycrystalline bulk ferromagnets^ the sources of 
spin disorder are related to lattice imperfections, e.g., 
point defects, dislocations, grain- and phase boundaries, 
or pores. These microstructural defects are accompanied 
by spatial variations of the materials parameters, for in¬ 
stance, the magnitude of the local saturation magnetiza¬ 
tion, exchange constant, or variations in the magnitude 
and/or direction of the magnetic anisotropy field. As a 
result, at a given value of the applied magnetic field, these 
features give rise to a deviation of the magnetization from 
the perfectly aligned state, in other words, they result in 
spin misalignment. On the other hand, inhomogeneous 
spin states in the bulk of a material (with V • M 0) go 
along with a magnetodipolar interaction field, which has 
an important impact on magnetic properties. 

The dipole-dipole interaction—one of the most funda¬ 
mental interactions in condensed-matter physics—is still 
the subject of current research. For instance, it gives rise 
to anomalous features in the ground-state correlations 
and in the spin-wave excitation spectrum of 2D spin sys¬ 
tems consisting of cold polar molecules;^ and it is vital 
for the understanding of spin-ice physics, where frustra¬ 
tion, dipolar ferromagnetic coupling, exotic field-induced 
phase transitions, and unusual glassiness are of relevance 
(see, e.g., Refs.an di for recent reviews). Recent exper¬ 
iments on PdFe islands on are square lattice using pho¬ 
toemission electron microscopy^ even indicate that pole 
interactions of higher order than the dipolar one are re¬ 
quired in order to understand the ground-state ordering 
features of such a system. 

While classical “standard” magnetometry provides 
only integral information about the magnetic state of 
the sample, scattering techniques, in particular, magnetic 
neutron scattering yield spatially and time-resolved in¬ 
formation about magnetic media. A further important 
difference between magnetometry and magnetic neutron 
scattering relates to the impact of the magnetodipolar in¬ 
teraction: the quantity of interest in an elastic magnetic 


neutron scattering experiment, the differential scattering 
cross section dTiM/dUl, depends in a twofold manner on 
the magnetodipolar interaction. First, and different from 
magnetization measurements, the interaction of the mag¬ 
netic moment of the neutron with the sample’s magneti¬ 
zation results in dipolar selection rules which are embod¬ 
ied, e.g., by the appearance of trigonometric functions in 
dYjM/dUl (via the Halpern-Johnson vector)i^ Second, the 
dipole-dipole interaction between the magnetic moments 
in the sample has a direct impact on its magnetization 
structure and therefore on the Fourier components of the 
magnetization. The former determine the magnetization 
(as measured by magnetometry) and the latter the prop¬ 
erties of dTiM/dUl. 

Traditionally, the influence of lattice defects on the 
magnetization of bulk magnetic materials is studied 
by analyzing magnetization curves in the approach-to- 
saturation regime. Early investigations by Brown,— us¬ 
ing the continuum theory of micromagnetics, provide an¬ 
alytic expressions describing different types of defects. 
From the neutron-scattering point of view, it is well 
known that the spin perturbations that are related to 
imperfections give rise to a strong magnetic scatter¬ 
ing signal—the so-called spin-misalignment scattering— 
along the forward direction (at small scattering angles) 4 
The size of perturbed regions is characterized by a field- 
dependent correlation length, which varies between about 
a few nanometers up to ~ 100 nm. It is the purpose of 
this paper to study the role of the magnetodipolar in¬ 
teraction on (real-space) magnetic correlations. In par¬ 
ticular, we aim to disentangle the twofold impact of the 
anisotropic dipole-dipole interaction on magnetic corre¬ 
lation functions obtained from neutron data. 

In our analysis we employ our previously developed mi¬ 
cromagnetic simulation methodology, which has provided 
fundamental insights into the magnetic small-angle neu¬ 
tron scattering (SANS) of various materials)^ The par¬ 
ticular strength of this micromagnetic approach is that 
it takes into account site-dependent magnetic interac¬ 
tions (for exchange, magnetic anisotropy, and saturation 
magnetization). This feature implies that the magnetic 
microstructure of a wide range of polycrystalline mag¬ 
netic materials such as single-phase nanocrystalline mag- 
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nets, magnetic nanocomposites, or magnetic particles in 
a nonmagnetic matrix can be studied. As a prime candi¬ 
date for a system with strong internal dipolar fields and 
nontrivial magnetization correlations, we have chosen for 
the present study porous ferromagnets (iron and cobalt), 
since here local variations in the saturation magnetiza¬ 
tion are at maximum. 


II. MICROMAGNETIC SIMULATION 
METHODOLOGY 

Our micromagnetic algorithm was originally developed 
for the computation of the magnetization distribution of 
magnetic nanocomposites and of the related magnetic- 
field-dependent SANS cross sectionsj^ The four standard 
contributions to the total magnetic energy (external field, 
magnetic anisotropy, exchange and dipolar interaction) 
are taken into account. In the simulations presented 
here, the sample volume V = 0.2 x 0.75 x 0.75 ^m^ 
was divided into N ^ 5 x 10® mesh elements, com¬ 
prising both pores and nanocrystallites. For the later 
comparison to experimental neutron data on nanocrys¬ 
talline inert-gas-condensed porous iron (Ref. 0) the vol¬ 
ume fraction of pores was chosen as P = 32 %, with 
randomly placed pore centers. Due to the flexibility of 
the mesh-generation algorithm, the shape of the pores 
can be controlled and was taken to be polyhedron-like. 
The pore-size distribution was assumed to be lognormalii 
with a median of 15 nm and a variance of 1.16, which 
yields a maximum of the distribution at 12 nm. The lo¬ 
cal saturation magnetization of each (iron) nanocrystal¬ 
lite was taken ^qMs = 2.2 T, which in conjunction with 
the above mentioned porosity value yields fioMs = 1.5 T 
for the entire sample. For the exchange-stiffness con¬ 
stant and the first cubic anisotropy constant of iron, we 
have, respectively, assumed values of A = 25pJ/m and 
ATi = 47kJ/m® (Ref. IT^ . The direction of anisotropy 
axes varies randomly from crystallite to crystallite. The 
energy-minimization procedure provides (at some partic¬ 
ular value of the applied magnetic field) the magneti¬ 
zation vector field M(r) = [Mx{r),My{r),M 2 {r)] of the 
sample on an irregular lattice. This distribution is then 
mapped onto a regular lattice, which permits us to cal¬ 
culate the magnetization Fourier coefficients and the en¬ 
suing neutron scattering cross section using Fast Fourier 
transformation. Further details can be found in Refs. 


III. RESULTS AND DISCUSSION 

Figure [1] displays the numerically computed spin mis¬ 
alignment Mx(r) = {Mx,My,0) in the vicinity of a 
nanopore together with the transversal Fourier coeffi¬ 
cients \Mx{q)\'^ and |Afy(q)p; these Fourier components 
will be used later on when discussing magnetic neutron 
scattering. The discontinuity (jump) of the magnetiza¬ 
tion at the pore-matrix interphase (/xq AM = 2.2 T) gives 



FIG. 1. (Color online) (left) Computed spin misalignment (at 
goHo = 0.6 T) around a pore {D = 12 nm) in a ferromagnetic 
iron matrix (2D cut out of a 3D simulation). Shown is the 
magnetization component M^(r) perpendicular to Ho || ez', 
thickness of arrows is proportional to the magnitude of Mx. 
Solid grey lines: magnetodipolar field distribution, (right) 
Corresponding magnetization Fourier components |Ma;(q)|^ 
and |My(q)p projected into the plane = 0. Bright colors 
correspond to “high” values and dark colors to “low” values 
of the Fourier components. Pixels in the corners of the images 
have q = 0.4nm“^. Logarithmic color scale is used. 


rise to a strong magnetostatic field and to the character¬ 
istic dipole-field-type spin texture. The symmetry of the 
Mx (r) distribution corresponds to the field of a magnetic 
dipole located at the center of the pore and aligned oppo¬ 
site to the external field direction, as depicted by the solid 
lines in Fig. [TJ Both transversal Fourier components are 
highly anisotropic, whereas the longitudinal one [M^p is 
isotropic (data not shown). We would like to emphasize 
that the angular anisotropy of both |Ma;p and |Mj,p is a 
consequence of the internal magnetodipolar field, which 
is the only long-range (nonlocal) anisotropic interaction 
in our micromagnetic modeling; neglecting this energy 
term results in |M 2 ;p and |Myp being isotropic^ 

However, in contrast to our previous results;^ where 
|Ma;p was found to be isotropic., the present simulations 
reveal an anisotropic |Ma;p Fourier coefhcient (which is 
enhanced along the held direction). This hnding is due 
to a corrected averaging procedure, which takes into ac¬ 
count magnetic fluctuations also along the e^, direction, 
which coincides with the direction of the incident neutron 
beam in a SANS experiment (see Fig. [T]). Specihcally, 
the computed (mapped on a regular grid) spin struc¬ 
ture of the sample, M.{x,y,z), is divided into thin slices 
i = l,...,Nx (with a typical thickness of 2.5nm). This 
results in a set of magnetization distributions 
and Fourier transformation then yields the g^). 
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The squared Fourier coefficients etc. are summed 

up and divided by the number of slices in order to obtain 
the averaged quantities. 

As mentioned above, the possibility to compute the in¬ 
dividual Fourier components allows one to compare the 
simulation results with experimental data for the spin- 
misalignment SANS cross section, which reads (for the 
scattering geometry where Hq || is applied perpendic¬ 
ular to the incoming neutron beam)j^ 

^^(q) = + \My\^cos^9 



-{MyU; + M*M^) sin 9 cos 9^ , (1) 

where V is the scattering volume, bn = 2.9x 10®A“^m“^, 
and 9 denotes the angle between the momentum-transfer 
vector q = {0,qy,qz) and (a* is the complex conju¬ 
gate of a). Note that dT,M/dfl represents the part of 
the total unpolarized SANS cross section dS/dVl which 
is exclusively due to transversal spin misalignment, with 
corresponding Fourier amplitudes and My; in other 
words, dTiM/dfl is obtained by subtracting the nuclear 
and magnetic SANS at complete saturation, dTi/dfl = 
^(|A^P +b'jj\Mz\'^ sin^ 9), from the measured dE/dfl at 

lower fields. __ 

The trigonometric functions and the MyM^-containing 
cross term in Eq. © are due to the neutron-magnetic 
interaction, while the Fourier components M^^y^z may 
additionally depend on the internal magnetodipolar in¬ 
teraction, e.g., the anisotropy of the dipolar interaction 
is embodied in the dependence of the M^^y on the angle 
9 (see Fig. [1]). Due to the complexity of this expression, 
the Fourier transform of dTiM/dfl, which one may call 
the correlation function of the spin-misalignment SANS 
cross section, does not represent, of course, the autocor¬ 
relation function of the magnetization. This is in contrast 
to the well-known result from nuclear scattering, where 
the nuclear scattering cross section 

^^(q) ^ j Cw(r) exp(-iqr) d^r (2) 

is equal to the Fourier transform of the autocorrelation 
function of the nuclear densityi^ 

C'Ar(i’) J AN{x) AN{x + r) d^x, (3) 

where AN{r) = N{r) — (N) denotes the so-called ex¬ 
cess scattering-length density, and (N) is the (constant) 
average scattering-length density, which only gives a con¬ 
tribution to dl^N/dft at q = 0. 

In order to compare the results for the correlation func¬ 
tion of the spin-misalignment SANS cross section, 

C'M(y,z)~y ^^^{qy,qz) exp{iqr)d‘^q, (4) 

with the autocorrelation function of the magnetization 
spin misalignment, C'sM(r), which is not decorated by 


FIG. 2. (Color online) Semiquantitative comparison between 
experimental datapi^ and micromagnetic simulations for the 
spin-misalignment SANS cross section of nanoporous iron 
(porosity: P = 32%; goHo = 0.6T). The respective scatter¬ 
ing signal at a saturating field of 1.83 T has been subtracted 
in both data sets. Logarithmic color scale is used. 


the neutron-magnetic interaction, we define the latter as 
follows 

CsM (r) J (x) (x + r) (5) 

Using the convolution theorem, this expression can be 
rewritten as 

C'sM(r) -- J (^\Mx{q)\'^ + |My(q)p) exp(fqr) d^q.{6) 

Note the additional difference between both correlation 
functions: Cm involves only a 2D integration over q 
space, because dT,M/dfl is experimentally accessible only 
in the plane perpendicular to the incoming neutron beam 
[small-angle approximation: q = whereas 

CsM is naturally obtained by a 3D Fourier transforma¬ 
tion [due to Mx^y = Mx,y{qx, <ly, <lz)]- Therefore, in order 
to compare both correlation functions, we have to em¬ 
ploy for the computation of Csm tbe same averaged 2D 
Fourier components |Ma;|^ and |My|^ that enter d'SM/dfl, 
i.e., the integration in Eq. is also only carried out over 
the y — z plane. 

Before discussing the correlations functions, we show 
in Eig. [5] a comparison between the experimental (spin- 
misalignment) SANS cross section of nanoporous ironic 
and the results of the micromagnetic simulations. The 
simulation data exhibit an overall good agreement with 
the experimental data and the anisotropic character of 
dT^M/dO, is clearly reproduced; this supports the validity 
of our micromagnetic simulation methodology. 

The correlation functions at 0.6 T, obtained using 
Eqs. ([ 3 ]) and ([ 5 ]), are depicted in Fig. [ 3 ] along the hor¬ 
izontal {z) and vertical (y) direction. One recognizes the 
existence of anisotropic correlations already for the auto¬ 
correlation function of the spin misalignment (not influ¬ 
enced by the interaction between neutrons and magnetic 
moments), which may be expected due to the long-range 
and anisotropic nature of the magnetodipolar interaction 
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FIG. 3. (Color online) Comparison between the normalized 
autocorrelation function of the spin misalignment Csm (solid 
lines) and the normalized correlation function of the spin- 
misalignment SANS cross section Cm (dashed lines) along 
different directions in the y — z detector plane (/ro^fo = 
0.6 T). The right images show the corresponding combi¬ 
nation of Fourier components, projected into the plane of 
the 2D detector: (autocorrelation) |Ma;p -|- \My\'^; (SANS) 
-I- cos^ d — {MyM^ + MyM^) sin 0 cos 6^. Pixels in 

the corners of the images have q = 0.4 nm~^. Logarithmic 
color scale is used. 


(see also Fig. [T|); the difference between both directions 
is significant (in particular for r = 30 — 40 nm) with Csm 
along the vertical direction being exclusively positive def¬ 
inite, while Csm along the horizontal direction intersects 
the r-axis at r = 20 nm and possesses a global mini¬ 
mum at r = 30 nm. The existence of anticorrelations in 
Csm around the particular r-value is a manifestation of 
the typical magnetization distribution Mj^(r) around a 
pore (see Fig. [T|), which is due to the configuration of 
the magnetodipolar field in the vicinity of such an inclu¬ 
sion. Namely, the perpendicular magnetization compo¬ 
nent changes its sign along the direction of applied field at 
a distance comparable with the pore diameter. Of course, 
the zeros and global minima of the correlation functions 
are dependent on the applied magnetic field. We point 
out (see also explanations above) that the difference be¬ 
tween the correlation functions corresponding to both 
approaches—neutron scattering and micromagnetics—is 
also significant, in particular, for r = 20 — 50 nm. 

The situation becomes even more interesting if the 
magnetic system possesses a second symmetry break¬ 
ing (in addition to the symmetry breaking caused by 
the external magnetic field). It can lead to the forma¬ 
tion of nontrivial 3D correlations, i.e., the correlation 
lengths of the magnetization distribution can be signifi¬ 
cantly different along the three coordinate axes. One can 
observe this effect, e.g., by introducing an antisymmet¬ 



FIG. 4. (Color online) Isosurfaces of the 3D magnetization 
Fourier components \Mx{ci)f and |Mj,(q)p of porous cobalt 
at an external magnetic field of go Ho — 0.6 T applied in z- 
direction. The image below schematically shows the assumed 
distribution of uniaxial anisotropy axes in the x — z plane. 


ric exchange contribution to the total energy, also called 
Dzyaloshinskii-Moriya interaction, favoring the forma¬ 
tion of spin canting 

Another example for a magnetic system with nontrivial 
3D correlations is a polycrystalline magnet with a non¬ 
isotropic {e.g., in-plane) distribution of anisotropy axes 
of the individual grains (see Fig. |4]). Here, we have sim¬ 
ulated this situation for the case of nanoporous cobalt. 
The choice of this material was motivated by the strong 
uniaxial anisotropy [Ku = 400kJ/m^ (Ref. E)]. The 
random 2D distribution of anisotropy axes is in the x — z- 
plane. All structural parameters including the porosity 
value are the same as we used for the simulations on 
porous iron. The local saturation magnetization and the 
exchange-stiffness constant of cobalt were, respectively, 
taken as g^Mg = 1.76T and A = 28pJ/m (Ref. [H). 

The simulation results are displayed in Fig. |T] as iso¬ 
surfaces of the 3D magnetization Fourier components 
\Ma,{qx,qy,qz)? and \My{qy„ qy, q:z)\^. One can clearly 
see that both Fourier components are strongly influenced 
by the magnetodipolar interaction, which results in a 3D 
clover-leaf-shaped pattern that was already presented in 
Fig. [Has a projection into the plane qx = 0. But the most 
important observation is that \Mx\'^ and |Mj,p are quali¬ 
tatively different (for a direct qualitative comparison one 
should rotate, e.g., the \Mx\^ component around tak¬ 
ing into account the relation between the Fourier compo¬ 
nents and the coordinate system). The magnetodipolar 
“structure” in |Mj,p is more pronounced, since the cor¬ 
responding coordinate vector By is perpendicular to the 
uniaxial anisotropy axes in the sample. This result shows 
that the magnetic correlations in porous cobalt with an 
in-plane distribution of anisotropy axes have a nontrivial 
3D character. The components can not be fully matched 
by rotation around q^, in contrast to the case of a 3D 
random distribution of anisotropy axes (see the results 
on porous iron above). 
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IV. SUMMARY AND CONCLUSIONS 


Using micromagnetic simulations we have computed 
the magnetization distribution of nanoporous ferromag- 
nets (iron and cobalt) which exhibits strong dipolar cor¬ 
relations due to nanoscale spatial variations of the mag¬ 
nitude of the magnetization. The results were used in 
order to compute the correlation function of the spin- 
misalignment SANS cross section, which is affected in a 
twofold manner by the magnetodipolar interaction, and 
the autocorrelation function of the spin misalignment, 
which by contrast depends only on the magnetization dis¬ 


tribution. Our approach was validated by comparing the 
simulation results with experimental neutron data. As 
a consequence of the long-range and anisotropic charac¬ 
ter of the internal magnetodipolar field, we find strongly 
anisotropic magnetic correlations. Our results demon¬ 
strate the importance of the magnetodipolar interaction 
for understanding magnetic neutron scattering. 
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